Small scale energy cascade of the solar wind turbulence 
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ABSTRACT 

Magnetic fluctuations in the solar wind are distributed according to Kol- 
mogorov's power law f~^^^ below the ion cyclotron frequency fd- Above this 
frequency, the observed steeper power law is usually interpreted in two different 
ways: a dissipative range of the solar wind turbulence or another turbulent cas- 
cade, the nature of which is still an open question. Using the Cluster magnetic 
data we show that after the spectral break the intermittency increases toward 
higher frequencies, indicating the presence of non-linear interactions inherent to 
a new inertial range and not to the dissipative range. At the same time the 
level of compressible fluctuations raises. We show that the energy transfer rate 
and intermittency are sensitive to the level of compressibility of the magnetic 
fluctuations within the small scale inertial range. We conjecture that the time 
needed to establish this inertial range is shorter than the eddy-turnover time, 
and is related to dispersive effects. A simple phenomenological model, based on 
the compressible Hall MHD, predicts the magnetic spectrum ~ /c"^/^"*"^", which 
depends on the degree of plasma compression a. 

Subject headings: solar wind turbulence. Hall MHD 

1. Introduction 



Solar wind, which is highly turbulent, represents a unique opportunity to investigate tur- 

bulence in natural plasmas via in situ measurements (iTu and Marsch 19951 : iBruno and Carbone 20051 ) . 
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In non-magnetized fluids, where the energy injection scale is far from the dissipation one, 
the intermediate scales (inertial range) are described by the universal power law Kolmogorov 
spectrum with s = 5/3. This law depends neither on the energy injection nor on the 
energy dissipation processes. 

In the solar wind, and in the interplanetary space in general, the mean free path roughly 
corresponds to the Sun-Earth distance and the usual dissipation via collisions is negligible. 
At the same time, in a magnetized plasma there is a number of characteristic space and 
temporal scales. Investigating solar wind turbulence at these scales is then a challenging 
topic from the point of view of basic plasma physics. Here we will focus our discussion on 
the ion scales, namely the ion inertial length Aj = c/upi and the ion cyclotron frequency 
fci = eB/iTLi. At these scales the fluid-like approximation of plasma dynamics, the usual 
Magnetohydrodynamic (MHD) description, breaks down in favor of a more complex descrip- 
tion of plasma. 

Solar wind turbulent spectrum of magnetic field fluctuations follows a ~ f^^^^ power law 
below the ion cyclotron frequency fd- For / > t he spectrum steepens significantly, but i s 
still described by a power law /~'^, with s G (2 — 4) (ILeamon et al. 19981 : ISmith et al. 20061 ). 



In situ solar wind measurements provide time series data, i.e. information on frequency 
in Fourier space. How to get any information on wave vectors? For fluctuations with 
velocities much smaller than the plasma bulk velocity V , the Taylor hypothesis is valid: 
the observed variations on a time scale 5t correspond to variations on the spatial scale 
5r = V6t. Therefore there is a direct correspondence between / and k spectra. If the 
solar wind turbulence was a mixture of linear wave modes, Taylor hypothesis would be well 
verified only below the spectral break: the bulk velocity is superalfvenic (V > Va, where Va 
being the Alfven speed) and so the low frequency fluctuations can be considered as frozen 
in plasma. However, whistler waves (with / > and phase speed > Va) do not satisfy 
the Taylor hypothesis assumption. In this study we assume that in the solar wind there are 
no whistler waves above the spectral break frequency. T his last assumption is supported by 
results recently obtained in the Earth's magnetosheath (IMangeney et al. 20061 ). 



Using the Taylor hypothesis the observed solar wind spectrum below the break is usu- 
ally attributed to the Kolmogorov's spectrum ~ k"^/^. Above the spectral break, the spec- 
tral steepening, ~ can be interpreted in two different ways. Some authors associate 



it to the dissipation range (ILeamon et al. 19981 : iLeamon et al. 19991 : iLeamon et al. 2000 



Bale et al. 20051 : ISmith et al. 20061). Others suggest that after the spectral break another 



turbulent cascade takes plac e (IBiskamp et al. 19961 : iGhosh et al. 1996 
Li et al. 20ol : ICaltier 20061 ). 



Stawicki et al. 2001 
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In ordin ary fluid flows the dissipation range is described by an exponential function 
( iFrisch 19951 ). While in the solar wind, a well-defined power law is observed after the break 
point and not an exponential. Note that power spectra, i.e. the second order statistics, 
completely describe Gaussian, or statistically independent, fluctuations. However, as is well 
known, fluctuations cannot b e described by a Gaussian statistics in the low-frequency part 
of the solar wi nd turbulenc e (IBruno and Carbone 20051 ). Deviations from Gaussianity, i.e. 
intermittency (IFrisch 19951 ). may be quantified by the flatness, the forth-order moment of 
fluctuations. 

In this paper we investigate the nature of magnetic fluctuations in the high frequency 
range of the solar wind turbulence, that is usually called dissipation range. We find that in 
this range the flatness increase with frequency. This is similar to what is going on in the 
low-frequency range. The presence of intermittency together with the well defined power law 
in the high-frequency part of the spectrum suggests another turbulent cascade rather than a 
dissipation range. This small scale cascade is observed to be much more compressible than 
the Kolmo gorov-like inertial ran ge, in agreement with the previous observations by the Wind 
spacecraft (iLeamon et al. 19981 ). We show that the energy transfer rate and intermittency 
are sensitive to the level of compressibility of the turbulent fluctuations in this range. Finally, 
we propose a simple phenomeno logical model, based on the compressible Hall MHD, which 
allows to explain the observed range of the spectral indices in the high frequency part of the 
solar wind spectrum. 



2. Turbulent spectra and intermittency 



In the present study we analyze the Cluster magnetic field data up to 12.5 Hz (or 
100/ci). We use 57 minutes of data during the interval 22:35:00-23:32:00 UT on 5 April 
2001, when Cluster was in the solar wind not connected to the bow-shock. During this 
period the interplanetary magnetic filed is B 7 nT, the ion density is rij ~ 3 cm~^, the 
plasma bulk velocity is V ~ 540 km/s, the ion temperature is Tj 33 eV, the Alfven 
speed is Va = B / y/JI^ff^rTi ~ 90 km/s, the ion skin depth Aj is about 130 km, the ion 
Larmor radius is pi ~ 80 km and the ion plasma beta is Pi = 2fioPi/B'^ ~ 0.8, pi being 
the pressure of the ions. All the data are from Cluster sp acecraft one, B is determined 
from the fluxgate magnetometer fFGM ) (IBalogh et al. 20011) . the other plasma parameters 
are provided by CIS/HIA instrument (iReme et al. 200ll ). Magnetic field fluctuations are 
measured by the search coils (SC) of the STAFF experiment with 0.04 s time resolution 
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( Cornilleau-Wehrlin et al. 2003 ). This instrument0 operates in the frequency domain (0.35 — 



12.5) Hz. We use FGM data to resolve the frequencies below 0.35 Hz. Results from Cluster 
have been compared with data from the Helios 2 satellite within the range (0.001 — l)fci- 
These data are selected within a high-speed stream V ~ 600 km/s when the satellite orbited 



at 0.96 AU flBruno et al. 20031 ). 



In order to analyze the magnetic field fluctuations in the solar wind we use the Morlet 
wavelet transform, defined as 

N-l 

m{r,t) = Y,B,{t,mt,-t)/r] (1) 

j=0 

which represents the transform of the i-th component of the magnetic field Bi{tj), a data 
time series with equal time spacing St and j = 0,...,N — 1. Here r is a time scale (the 
corresponding frequency is / = l/r), while ip{u) = 2-'^/^7r^-'^/'' cos(ti;oM) exp[— 7x^/2] is the 
Morlet wavelet, where uq = 6. The square of a wavelet coefficient gives a "quantum of 
energy" of magnetic fluctuations on a scale r at a time t. Thus we define the power spectral 
density of Bi as 

9f)f ^"^ 

5anTVHz] = — 5^|W.(r,t,)r. (2) 

j=0 

Fig.[T]shows the total power spectrum density of the magnetic fiuctuations S = J2i=x y z^i 
(solid lines) and the spectrum of magnetic field modulus fiuctuations 5|B| (dashed-dotted 
lines), which corresponds to the spectrum of parallel fiuctuations of magnetic field S\\. In 
agreement with previous observations, the solar wind spectra follow well defined power laws 
with a break around ~ 0.3 Hz, in the vicinity of — 0.1 Hz. Taking into account 
the Taylor hypothesis, fb corresponds to a spatial scale V/ fb — 1800 km ~ 15Ai and to 
the normalized angular wavenumber kXi ~ 0.4 {kpi ~ 0.3). Below the break, the Helios 
spectrum S ~ f^^'^^ is fitted roughly by Kolmogorov's law. The Cluster FGM spectrum 
covers a (0.02 — 0.5) Hz frequency range, the lower resolved frequency here is determined 
by the length of the consi dered time period, 57 minut es, and by a cone of infiuence of the 



Morlet wavelet transform (ITorrence and Compo 19981 ). Even if the frequency range is not 



that large, this spectrum can be described by a power law 5 ~ / Above the break, the 
spectrum follows a S ^ f~^'^ law: t he spectral index 2.6 lies in the usual observed range. 



between 2 and 4 (jSmith et al. 20061 ). As regards the spectrum (Sy, it is remarkable, that 



^STAFF-SC measurs magnetic fluctuations above 0.1 Hz. As the Cluster spin frequency is 0.25 Hz, there 
are artefacts at 0.15 and 0.35 Hz in the solar wind. We note that in the magnetosheath, where the level of 
magnetic fluctuations is much higher than in the solar wind, this problem is almost unnoticeable. 
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after the spectral break, the compressible magnetic fluctuations become much more impor- 
tant than in the low frequency part of the spectrum. This can be quantified by the ratio 
(S), that is 0.02 — 0.05 in the law frequency part and 0.26 in the high frequency part of 
the spectrum. This means that nearly 30% of the magnetic fluctuations above the spectral 
break are compressible. 

Let us now consider statistical properties of turbulent fluctuations in the high-frequency 
part of the spectrum. Fig. [2] shows the probability distribution functions (PDFs) of the nor- 
malized variables = yVx{T,t)/crx, with the standard deviation = (Wa;(r, t)^)'^''^, at 
three different time scales r; the dashed lines indicate the corresponding Gaussian fits. Ap- 
parently, the smaller the scale, the more the PDFs deviate from Gaussian. The PDFs of the 
other components of the magnetic field present a sim ilar behaviour. The dependence of nor- 



malized PDFs on scale is a signature of intermittency (IFrisch 19951 : iBruno and Carbone 20051 ) 



Intermittency may be quantified by the flatness. Fig. [3] shows the spectrum of flatness 
F{f) = / {W^y of fluctuations. We can see that in the low-frequency part F 

increases from the large-scale Gau ssian value F = 3. This means that the usual Alfvenic 



cascade (IBruno and Carbone 20051 ) proceeds towards smaller scales building up phase corre- 
lations among fluctuations. In the vicinity of fd, Helios indicates an increase of the flatness 
and Cluster/FGM shows a slight decrease. This disagreement can be due to (i) aliasing at 
the highest frequencies of Helios; (ii) artifacts on FGM at frequencies around the satellite 
spin frequency 0.25 Hz; or to some other physical reasons which are outside of the scope of 
the present paper. Here we will focus our discussion on the higher frequencies, resolved by 
Cluster/STAFF-SC, where the flatness again displays a power-law dependence on /. This 
scaling behavior is inherent to non-linear dynamics and refers to a new kind of "inertial 



region 

In non-magnetized fluid, a rapid exponential increase of the flatness is observed in the 
near-dissipation rang e since only the strong est fluctuations survive while the others are de- 



stroyed by viscosity (IChevillard et al. 20051 ). After such increase the flatness saturates. In 
our case, the saturation is not observed and there is a power law increase of flatness. At 
the same time strong coherent structures which give an important increase of intermittency 
represent only (2 — 6)% of the fluctuations at each scale (see the insert of Fig. [3]). This esti- 
mation has been made using a method proposed by Farge (1992), named local intermittency 
measure analysis. It consists in separating the intermittent events which participate to the 
non-Gaussian tails of the PDF's from the Gaussian fluctuations. At each scale r we chose 
a threshold ( for the energy of magnetic fluctuations and we select the wavelet coefficients 
which verify |>V(r, t)p < (. With the new set of coefficients we recalculate the flatness; if 
it is larger than 3 we decrease the threshold until the flatness reaches the Gaussian value. 
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At the end of the procedure, the fraction of the non-Gaussian coefficients gives an estimate 
of the percentage of the intermittent structures at each scale. Note, that this analysis only 
gives the upper limit of the structures because one intermittent event is described normally 
by more than one wavelet coefficient. 



3. Role of compressible fluctuations 

The results described in the previous section suggest that above the spectral break 
frequency fb there is a non-linear compressible cascade and not a dissipation range. In this 
section we show that the level of the plasma compressibility in this range can be at the origin 
of a n on- universality of statistics of the magnetic fluctuations. To show this we consider two 
7 min time periods on April 5, 2001, with diflerent ion plasma (3i. 

Fig. m^a) shows the total power spectrum density of the magnetic fluctuations within 
the high frequency range S (solid line), and the spectrum of parallel fluctuations of magnetic 
fleld S\\ (dashed-dotted line) for a time period which starts at 23:03:02 UT, when j3i ~ 0.5. 
Fig. m^b) shows S and S\\ for a period starting at 22:35:02 UT, when (3i ^ 1.5. The straight 
lines indicate power law flts: S ~ f""^'^^ for the flrst time period and S ~ Z"^'^" for the 
second. Fig. Hl^c) shows the ratio S\\/S for the two considered time periods. One can see 
that during the flrst time period (with jSi ~ 0.5) the energy of compressible fluctuations is 
less than 20% of the total energy of the fluctuations (see the solid line), while for the second, 
it is around 30% (see the dashed line). Fig. Hl^d) shows the flatness F{f) calculated within 
the two time periods (here again the solid line corresponds to the period with smaller /3j and 
dashed line, to the period with higher Pi). 

These observations indicate that for higher the level of compressibility increases, 
leading to a steepening of the magnetic power law and an increase of intermittency. To 
conflrm these results more intervals for wider plasma parameters should be analyzed. 



4. Nature of the small scale cascade 

How can we explain the small scale energy cascade in the solar wind turbulence, sen- 
sitiv e to the compressibility o f the fluctuations? There are several interpretations. Follow- 



ing (ISchekochihin et al. 20071 ). at spatial scales between the ion and electron characteristic 
scales, there is a superposition of (undamped) Kinetic Alfven Wave (KAW) and entropy cas- 
cades. KAW waves are expected to have low frequencies / <C fd and nearly perpendicular 
wave vectors with respect to the mean fleld, k± S> k\\. 



- 7- 



Another possible explanation of the small scale cascade is a superposition of non-linear 
fluctuations (not waves), which exchange energy on times smaller than the usual eddy- 
turnover (or non-linear) time. At spatial scales at the order of Aj and at frequencies of the 
order of fd one should take into account the Hall effect in Faraday's equation 



dB „ 



m ■ 

VxB -J X B + r^J 



(3) 



J = cV X B/47r being the current density and t] a scalar dissipation coefficient. This equation 
contains three different physical processes, with three associated characteristic times. By 
introducing for each length-scale i the mean value of density pe, magnetic intensity 
and velocity ue, we can define an eddy-turnover time T/vl ~ ^/ue, a characteristic time 
associated to the Hall effect ~ Pii'^/Bi and a dissipative time Td ~ i'^/rj. At large scales 
the main effect is played by the term V x (V x B), which describes the low- frequency Alfvenic 



cascade (IBruno and Carbone 20051 ) realized in a time T^l- We conjecture that approaching 
the scale of the spectral break, Th appears to be of the order of T/vl and beyond the break, 
the nonlinear energy transfer is realized in a time ~ < Tnl,T£,. 

In a compressible fluid the energy balance equation must be expressed in terms of 
energy densities (i.e., energy per unit volume) and not in terms of specific en ergy (i.e., per 



unit mass), in order to take into account density fluctuations (see (IFleck 19961 ) for example). 



The mean volume rate of energy transfer on Hall times in a compressible fluid is therefore 

~ -TT^- (4) 

Assuming equipartition between kinetic and magnetic energies peuj ~ -B|, the energy- 
transfer rate results to be proportional to 

ey^^r^-^ (5) 

Following von Weizsacker (1951), the ratio of the mass density p at two successive levels 
of the hierarchy is related to the corresponding scale size I by the following equation 



-3a 



(6) 



where |a| is a measure of the degree of compression at each level v (larger v meaning 
larger length-scale), and ranges from a = for no compression up to |a| = 1 for isotropic 
compression (3|a| is a number of dimensions in which the compression takes place). So, 
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using this density scaling, pe ^ i and assuming a constant spectrum energy transfer rate 
we have 

Be ~ p'/^'/' ~ (7) 
Therefore, the spectral energy function is 

E(k) ~ 4^ ~ fc-^/3+2a_ 

k 

In the incompressible plasma limit (a = 0), this phenomenology predicts a ~ 
spectrum. Such a spectrum has been observed b oth in direct numeric al simulations of an 



incompressible Electron MHD turbulent system (IBiskamp et al. 19961) and in the EMHD 
limit of the incompressible Hall MHD shell model (jcaltier and Buchlin 2007 ). 



In the case of isotropic compression toward smaller scales {a = 1) which can take place 
in the interstellar medium, the spectrum is E(k) ~ k~^^^. If the isotropic compression is 
going on toward larger scales (a = —1) the spectrum will be E(k) ~ k~^^^^. In the case, 
analyzed here, the power law spectrum of magnetic fluctuations follows the ~ A;~^'^-law, 
that corresponds to a ~ —0.14, i.e., the scaling relation for the density is pe ~ 

This simple phenomenological model allows to explain the variations of the spectral 
index of the high frequency part of the solar wind spectrum from —4 to —2 by different 
degree of plasma compression with a between —5/6 and 1/6. 

However, the described model have to be improved in order to take into account (i) 
the space anisotropy {k± ^ k\\) that appears in a plasma with a mean field and (ii) possible 
differe nt scaling for velocity and magnetic field which can be found in the compressible Hall 



MHD (jServidio et al. 20071 ). Besides, the interpretation of the small scale energy cascade 
in the Hall MHD frame is supported by the observations of a clear correlation between the 
spectral break frequency fb and the Doppler-shifted wavevector of the ion inertia length /a, 
(Leamon et al, 2000). 

Note that the solar wind spectral break frequency fb is usually observed to be lower 



than (iLeamon et al. 19981 ). It can be explained by the fact that the Hall effect starts 
acting at scales larger than Aj. In our case, for example, fb/ f\i — 0.1 and the scale of the 
spectral break is ~ lOAj. 

Introducing the Hall Reynolds number as the ratio of the non-linear term over the Hall 

term 

nn — (9) 

cm,- Bf 
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one can see that Rh — 1 (i.e. Tnl — Th) at a scale £, that is not necessarily equal to Aj, 

— =47r^ 

here and are the mean values of the Alfven speed and the ion inertial length at scale 
i respectively. 



5. Conclusion 

In this paper we investigate small scales turbulent fluctuations in the solar wind (i.e. at 
frequencies above the spectral break at the vicinity of fd)- Taking into account the Taylor 
hypothesis {£ = V/ /) the frequency domain above the spectral break covers ~ (40 — 2000) km 
space scales that corresponds to (0.3 — 15)A.j. In this range we found evidences for strong 
departure from Gaussian statistics and the presence of intermittency while the spectrum 
presents a well-defined power law. Both the presence of a power-law spectrum and the 
absence of global self-similarity, seems to be quite in contrast with the role of "dissipative 
range" . 



In usual fiuid turbulence, the dissipative range (IFrisch 19951 ) starts with a rough expo- 



nential cutoff; in the near dissipation range the intermitten cy increases as far as th e Gaussian 



fiuctuations dissipate faster than the coherent structures (IChevillard et al. 20051 ): then the 
fiuctuations become self-similar, the singularities being smoothed by dissipation. In the solar 
wind turbulence we observe a completely different picture. After the spectral break in the 
vicinity of fd the flatness increases as a power law indicating that non-linear interactions 
are at work to build up a new inertial range. 

This small scale cascade is much more compressible than the lower frequency Alfvenic 
cascade. This sudden change in nature of the turbulent fluctuations can happen due to 
a partial dissipation of magnetic fluctuations at the spectral break: th e left-hand Alfveni c 
fluctuations with k\\ ^ k± are damped by the ion-cyclotron damping (iGhosh et al. 19961 ). 



Above the break, however, a new magnetosonic cascade takes place up to the electron char- 
acteristic scales. This energy cascade is seems to be dominated by the fluctuations with 



> (ISahraoui et al. 2009 : Mangenev et al. 20061 ). 



We found, as well, that the plasma compressibility controls the statistics of the magnetic 
fleld fluctuations. Preliminary results show that the increase of the level of the plasma 
compressibility leads to the spectrum steepening and increase of the intermittency. 



To explain this small scale compressible cascade we introduce here for the flrst time 
a simple phenomenological model based on the compressible Hall MHD: we flnd that the 
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magnetic energy spectrum follows a E{k) ~ j^-^/^+'^a j^^^^ depending on the degree of plasma 
compressibility a. While far from a complete model of small scale turbulence, this simple 
model can explain observed variations of the spectral index within the high fr equency part 



of the solar wind turbulent spectrum (ILeamon et al. 19981 : ISmith et al. 20061 ) by different 
degree of plasma compressibility. 

We acknowledge discussions with A. Mangeney, R. Grappin, R. Bruno and S. Servidio. 
We thank N. Cornilleau, E. Lucek and I. Dandouras for providing Cluster data. OA ac- 
knowledges the support of CNES and European network TOSTISP. 
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Fig. 1. — The total magnetic field power spectral density <S and the spectrum of compressible 
magnetic fluctuations S\\ (dashed-dotted line) measured by Helios 2 (up to 0.08 Hz) and by 
Cluster (up to 12.5 Hz). Straight lines refer to power law fits. Vertical dotted lines indicate 
the ion cyclotron frequency fd and the Doppler-shifted ion inertial length /;^.. 
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Fig. 2. — The probability distribution functions (PDFs) of the normahzed variable hxij) at 
three time scales: r = 2 s (upper panel), r = 0.7 s (middle panel), and r = 0.2 s (bottom 
panel). For each panel a Gaussian function is plotted for reference as a dashed curve. 
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Fig. 3. — The flatness of the magnetic field component calculated using Cluster data (solid 
line) and Helios 2 data (dashed-dotted line); the horizontal solid line indicates the Gaussian 
value F = 3. The insert gives the percentage of intermittent events which participate to the 
flatness in the high frequency part of the spectrum. 
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Fig. 4. — Properties of magnetic fluctuations above the spectral break for different plasma 
Pi] (a) S and ^Sy for a 7 minutes time period with /5j ~ 0.5; (b) the same as (a) but for a 7 
minutes time period with /3j ~ 1.5; (c) level of compressible fluctuations for period (a), solid 
hue, and (b), dashed line; (d) flatness of the fluctuations of component calculated for the 
period (a), solid hue, and for the period (b), dashed line. 



